Oscillatory thermal instability — the Bhopal disaster 

and liquid bombs 



Rowcna Ball 

Mathematical Sciences Institute, The Australian National University 
Canberra ACT 0200 Australia 



Abstract 

Thermal runaway reactions were involved in the Bhopal disaster of 1984, in 
which methyl isocyanate was vented from a storage tank of the liquid, and oc- 
cur in liquid peroxide explosions, yet to date there have been few investigations 
into the mechanism of thermal runaway in such liquid thermoreactive systems. 
Consequently protocols for storing thermally unstable liquids and deactivating 
liquid bombs may be suboptimal. In this work the hydrolysis of methyl iso- 
cyanate and the thermal decomposition of triacetone triperoxide were simulated 
using a gradientless, continuous-flow reactor paradigm. This approximation en- 
abled stability analyses on the steady state solutions of the dynamical mass 
and enthalpy equations. The results indicate that thermal runaway in both sys- 
tems is due to the onset of a large amplitude, hard thermal oscillation initiated 
at a subcritical Hopf bifurcation. This type of thermal misbehaviour cannot 
be predicted using classical ignition theory, and may be typical of liquid ther- 
moreactive systems. The mechanism of oscillatory thermal instability on the 
nanoscale is elucidated. 

Keywords: Thermal runaway, Non-classical ignition, Bhopal disaster. Methyl 
isocyanate, Organic peroxide explosions 



1. Introduction 

The causes of the world's worst industrial disaster at Bhopal on 2-3 Decem- 
ber 1984 are still being debated in the international media, more than 25 years 
after the event. Was it caused by neglect, parsimony, or procrastination by 
Union Carbide on training, safety and maintenance? Corruption, sabotage and 
cover-up? Inadequate government regulation? Any or all of the above may well 
have helped set up the worst possible scenario — for it could not have been any 
worse — but they are contributing factors rather than causes. (A brief account 
of the disaster is given in the Appendix.) 
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The primary cause of the thermal runaway that led to the venting of a lethal 
mist of methyl isocyanate (MIC) over Bhopal city was physicochemical. In this 
work I present a stability analysis of a simple dynamical model for the MIC- 
H2O reacting system, revealing oscillatory thermal misbehaviour that cannot be 
predicted using classical ignition theory. A similar instability is shown to govern 
the explosive thermal decomposition of the organic hydroperoxide triacetone 
triperoxidc, in liquid solution. The provenance of oscillatory thermal instability 
on the nanoscale is elucidated, and shown to lie in the ability of the reactant 
molecules to store energy in the internal molecular motions — in other words, 
the heat capacity. 

Despite the enormity of the Bhopal disaster little or no research has been 
published that elucidates the fundamental physico-chemical cause of thermal 
runaway in liquid reactive systems such as MIC hydrolysis. In terms of achiev- 
ing the Millenium Development Goals (MDG) it seems rather important that 
thermoreactive processes in liquids are thoroughly investigated and the knowl- 
edge disseminated widely. Given the the horrific legacy of the disaster, the long 
term adverse health effects in children and adverse reproductive eflfects in women 
of MIC exposure that have been well-documented (Mishra et al., 2009), such 
knowledge is relevant to the MDG of Child Health and Maternal Health. More 
generally people have a right to expect that thermally unstable and hazardous 
liquids are stored safely. 

The Millennium Development Goals Report (www.im.org/millcnniumgoals, 
2010) highlights the challenges posed by conflicts and armed violence to human 
security and MDG achievements. A new and growing threat to people's security 
is the use of liquid peroxide explosives by terrorists. The ingredients for making 
such bombs are cheap and widely available and they cannot be detected by 
metal detectors and nitrogenous explosives detectors, or distinguished from hand 
lotion by x-ray machines. Liquid peroxide-based explosives were used in the 
suicide attacks on the London transit system in 2005, which killed 56 people, 
and the terrorists convicted of the foiled 2006 transatlantic aircraft conspiracy 
had plotted to blow up a number of planes using liquid peroxide explosives. 
(Many more accounts of peroxide misuse incidents are easily found on the web.) 
To this day there are severe restrictions on carrying liquids through security 
barriers at most airports. It seems grimly inevitable that the use of liquid 
peroxide explosives as mass murder weapons will increase. Knowledge of their 
fundamental mechanism of action may help to counter their use. 

The open literature on theoretical and experimental validations of thermal 
runaway criteria and parametric sensitivity in batch reactors and storage tanks 
was summarized by Velo et al. (1996). In defining critical conditions they, along 
with other authors cited therein, begin with the assumption that storage tanks 
can be modelled as well-stirred batch reactors with linear thermal coupling to 
the environment. However batch reactors have no non-trivial steady states, 
and there is no general theory for determining whether a thermal excursion will 
grow or decay. It is shown in this work that a simple model with nonequilibrium 
steady states that is also spatially homogeneous — the continuous-fiow stirred 
tank reactor (CSTR) paradigm — can provide great insight into thermoreactive 
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instabilities in liquid systems, and provide fundamental causative information 
that cannot easily be extracted from detailed numerical simulations that include 
convective motions. 

In section 2 I describe the chemical reactions and provide the relevant data, 
taken from the literature, for the physical properties of the reactants and ther- 
modynamic and kinetic parameters. The CSTR paradigm is described in section 
3 and the equations are given using dimensionally consistent units and in terms 
of dimensionless variables and parameter groups. In section 4 the results of 
numerical stability analyses of the equations are given, where numerical values 
of the parameters for MIC hydrolysis and for triacetone triperoxide thermal de- 
composition in solution were used in the equations. Some points regarding the 
applicability of the CSTR paradigm are discussed in section 5, and the nanoscale 
aspects of oscillatory thermal rimaway are elucidated through examining the 
behaviour in the limits of the two timescales of the relaxation oscillation. A 
summary of the conclusions is given in section 6. 

2. Chemistry and data 

2.1. MIC hydrolysis 

Isocyanates hydrolyse exothermically to the corresponding amine and carbon 
dioxide. In excess water isocyanates react exothermically with the hydrolysis 
product amine to form the disubstituted urea (Saunders and Slocombe, 1948; 
D'Silva et al, 1986). With MIC the product is N,N-dimethyl urea and the 
reaction sequence is as follows: 

CH3NC0(i) + H20(i) '-^^ CH3NH2(aq) + C02(aq) (Rl) 
CH3NC0(i) + CH3NH2(aq) ''^^ CH3NHCONHCH3(aq) . (R2) 

For reaction R2 Ai?2(298K) = — 174.6kJ/mol and for the sequence overall 
Ai7toi(298K) = -255kJ/mol (Lidc, 2008-2009). 

A chemical analysis of the residue in the MIC storage tank (Tank 610) at the 
Union Carbide plant at Bhopal, sampled seventeen days after the event, found 
a variety of MIC condensation products, mainly the cyclic trimer (D'Silva et al., 
1986). However, experiments indicated that these condensations must have been 
initiated at temperatures and pressures well above the normal boiling point of 
MIC, so for modelling the initial thermal runaway these reactions need not be 
considered. No kinetic data are available for reaction R2 so only reaction Rl 
is used in the model. It will be seen from the results that reaction Rl alone is 
sufficient to induce thermal runaway. Relevant physicochemical data are given 
in table 1. 

2.2. Thermal decomposition of triacetone triperoxide (TATP) 

Triacetone triperoxide, a cycle trimer, is an explosive made by mixing ace- 
tone and hydrogen peroxide, both of which substances are legal, cheap and 
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Table 1: Physical, kinetic, and thermochemical parameters for MIC hydrolysis. 



Aiolecular weight MIC 


57.051 






Specific heat capacity C°(298) MIC 


1959 J/ (kg K) Poling 


et al. 




(2008) 




Specific heat capacity C°(298) H2O 


4181 J/(kgK) 






BoiUng point MIC at 1 atm 


38.3°C 


Lide 
2009) 


(2008- 


Density MIC at 25° C 


0.9588 g/cm^ 


Lide 
2009) 


(2008- 


Rl reaction enthalpy 


80.4kJ/mol 


Lide 

2009)t 


(2008- 


Rl activation energy 


65.4kJ/mol 


Castro 


et al. 




(1985) 




Rl pseudo first order frequency factor 


4.13e08/s 


Castro 
(1985) 


et al. 



'''From standard enthalpies of formation at 298 K. 



readily available over the counter. Pure TATP is a white crystalline powder 

that is soluble in organic solvents. The thermal decomposition of TATP does 
not involve combustion; the main reaction products are acetone, some carbon 
dioxide, and ozone (Eyler et al., 2000; Oxley et al., 2002). Its high explosive 
power is in part due to the large entropy increase of the formation of four gas 
molecules from one condensed-phase molecule (Dubnikova et al., 2005). Rele- 
vant parameters for the thermal decomposition of TATP in toluene are given in 
table 2. 



3. The CSTR paradigm 

The spatially homogeneous flow reactor, or reacting mass or volume, in 
which a reactant undergoes a flrst order, exothermic conversion is a simple but 
elucidatory model for thermoreactive systems when it is appropriate to ignore 
convection, because as a dynamical system it has non-trivial steady states that 
can be analysed for stability. The dynamical mass and enthalpy equations may 
be written as 

=M^e-^/«^c + F{cf - c) (1) 

at 

MCr^ = {-AH)Mze-^/'^^c + F{CfTa - CrT) - L{T - T„). (2) 

The symbols and quantities are defined in table 3. For numerical and com- 
parative reasons it is more convenient to work with the following dimensionless 
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Table 2: Physical, kinetic, and thermochemical parameters for thermal decomposition of 

TATP in toluene. 

Molecular weight 

Specific heat capacity C° (298) toluene 

Boiling point toluene at 1 atm 
Density of toluene at 298 K 
Reaction enthalpy 

Activation energy 

Frequency factor 

Feed concentration of TATP 
^Depending on reaction products. 

system corresponding to equations (1-2): 

xe-^'^+ef{^Ua-u)-l{u-Ua), (4) 

where the dimensionless groups are defined in table 3. Numerical analysis of 
equations (3-4) was carried out using values of the dimensionless groups calcu- 
lated from the data in tables 1 and table 2 and assigned values of the inverse 
residence time /, heat loss coefficient i, and inflow concentration c/. 

4. Results 

4.1. MIC hydrolysis 

From equations 3 and 4 in the steady state we can define the rate of reactive 
heat generation as the nonlinear term 

r,^/e-V7(e-V« + /), 
and the rate of non-reactive cooling as the linear term 

r^ = -u{ef + (.) + Ua{e^f + l). 

From classical ignition theory the reacting mixture self- heats if rg exceeds r;. 
Thermal runaway occurs if rg exceeds r; beyond a system-specific threshold; 
for the hydrolysis of MIC this is taken as the boiling point of MIC. These 
rates were computed using data from table 1 and plotted in figure 1, where 
the temperature is labeled in dimensional units. We see that the system self- 
heats until the reaction temperature T reaches the steady state temperature of 



222.2356 g/mol 

1698.25 J/(kg Haling et al. 
(2008) 

110.8 °C 

866.9 kg/m^ 

330- Dubnikova 
420kJ/molt et al. (2005) 
178.52 kJ/molEyler et al. 
(2000) 

9.57el6/s Eyler et al. 
(2000) 

2mol/kg 



da; 
du 
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Table 3: Quantities, definitions, and units. 



A 


reaction frequency 




c 


c{t), concentration of reactant 


mol/kg 


c 


inflow reactant concentration 


mol/kg 


Cf 


specific heat capacity of reaction mix- 


J/kgK 




ture 




Cf 


specific heat capacity of the inflow 


J/kgK 




stream 




E 


activation energy 


J/mol 


F 


fiow through rate 


kg/s 


AH 


reaction enthalpy 


J/mol 


M 


mass of reaction mixture 


kg 


R 


gas constant 


8.314 J/mol K 


t 


time 


s 


T 


T{t) reaction temperature 


K 


T 


ambient temperature 


K 


L 


heat loss coefficient 


W/K 


e 


CrE/cf{-AH)R T 


tA 


f 


F/MA u 


RT/E 


7 


Cf/Cr Ua 


RTa/E 


£ 


LE/cfMA{-AH)R x 





~305 K at which the heating and cooling rates are balanced. Since the boiling 
point of MIC is 312 K, according to this diagram the Bhopal disaster did not 
happen. On the basis of this diagram we would not expect a thermal runaway 
to develop, even when the ambient temperature is allowed to drift slowly up to 
292 K. 

However thermal balance diagrams such as that in figure 1 can be danger- 
ously misleading because they infer stability rather than assess stability rigor- 
ously, although such diagrams are often used in chemical reactor engineering. 
The steady states, periodic solutions, and stability analysis of equations (3 4) 
were computed numerically (Doedel, 2010) and yielded a dramatically different 
picture of the the thermal stability of MIC hydrolysis. Figure 2 shows a bi- 
furcation diagram in which the steady state temperature and the temperature 
amplitude envelope of periodic solutions are plotted as a function of Ta- 

The steady state is stable at « 286 K, the temperature at which the tank 
of MIC had been held for several months. As Ta is quasistatically increased 
the reaction temperature T increases slowly, but at = 290.15 if the stability 
analysis flags an abrupt change in the nature of the solutions. At this point 
the steady state solutions become unstable at a Hopf bifurcation and the hard 
onset of a high amplitude thermal oscillation ensues. Clearly, at Ta = 292 K we 
have catastrophic thermal runaway, contrary to the prediction given by figure 1. 
(In the resulting superheated fluid the exothermic condensation reactions would 
increase the temperature even further.) 
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280 290 300 310 320 K 



Figure 1: Rates of reactive heat generation Vg (red) and heat loss r; (blue) versus T from 
equations (3-4). Ua = 0.0379 (corresponding to Ta = 292 K), f = 1.7, £ = 700, e = 10. 

This is quite different from classical ignition of a thermoreactive system, 
which occurs at a steady-state turning point. The dynamics of oscillatory ther- 
mal runaway can be understood by studying the close-up of the region around 
the Hopf bifurcation Hi shown in figure 3. Hi is subcritical and the emergent 
limit cycle is unstable. The amplitude envelope of the unstable limit cycles is 
marked with a thin dotted line; they grow as Ta is decreased. At the turning 
point * of the periodic solution branch the limit cycles become stable. Thermal 
runaway may occur if there are significant perturbations while Ta is within the 
regime H«-i/i, and it must occur when Ta drifts above Hi. In principle the 
rapid thermal excursion takes the system to the amplitude maximum of the 
stable limit cycle. In reality the reactant and products have vaporised, the pres- 
sure has soared, the peak temperature is far above the boiling point of MIC, 
and the system must vent or explode. But it must be emphasized that the 
thermal runaway is due to oscillatory instability rather than classical ignition 
at a turning point. 

The presence of oscillatory instability is all-pervasive and dominant in this 
system. This can be appreciated by inspection of figure 4, a plot of the loci 
of the steady state turning points and the Hopf bifurcations of equations (3-4) 
over the two parameters Ua and the inverse residence time /. In the filled region 
thermal runaway will always be oscillatory. The bistable regime, indicated by 
the dashed line, occurs at very high flow rates (short residence times) . However, 
classical thermal runaway at a steady state turning point does not occur because 
the oscillatory instability is still present and dominant. 

Two computed time series for F — 0.0016 kg/s are compared in figure 5. 
In the upper plot Ta = 308.4 K and the oscillations decay to a stable steady 
state. However, the onset of thermal instability is violent: in the lower plot 
Ta = 308.5 K and the transient does not decay but swings into sustained high 
amplitude relaxation oscillations with a period of about 166 s. 

The behaviour of the system under a slow upwards drift of the ambient 
temperature can be simulated easily; a time series with dTa/dt — 0.02°C/s 
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Figure 2: Bifurcation diagram. Stable steady states are plotted with solid blue line, unstable 
steady states with dashed red line, and the amplitude envelope of periodic solutions is marked 
with thin dotted magenta line. Hi and H2 label the Hopf bifurcation and the large * marks 
the change in stability of the limit cycles, f = 1.7, £ = 700, e = 10. 



is shown in figure 6, which confirms the abrupt onset of the instabiUty. Of 
note is the decay in amphtude of the oscillations as the thermostat temperature 
increases; physically this occurs because the reactant is consumed faster than 
it is supplied as the temperature increases. 

4.2. Triacetone triperoxide 

The steady states, periodic solutions, and stability analysis of equations (3- 
4) were computed using the data for TATP thermal decomposition in table 2 
and results are shown in figures 7 and 8. 

In figure 7 the loci of the Hopf bifurcations are plotted over Ta and F, 
and a point within the oscillatory region was selected to compute the time 
series in figureS. The system exhibits violent relaxation oscillations, suggesting 
that explosive thermal decomposition of TATP is initiated at the onset of this 
oscillatory behaviour rather than by classical ignition. 
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Figure 6: Computed time series for MIC liydrolysis with drift in tliermostat temperature of 
0.02°C/s. F = 0.0016 kg/s, L = 560 W/K. 
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Figure 7: The locus of Hopf bifurcations for the triacetone triperoxide system. 
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Figure 8: Strong relaxation oscillations in the TATP thermal decomposition. The time series 
was computed for Ta = 460K, F = 5e-04kg/s. 
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5. Discussion 

The tendency for oscillatory thermal runaway may be typical of oxothcrmi- 
cally reactive organic liquids. In the work of Ball and Gray (1995) the hydration 
of 2,3-epoxy-l-propanol in a CSTR was found to exhibit similar non-classical 
thermal misbehaviour. Here it is shown that the thermal runaway that led to 
the Bhopal disaster may have been initiated at an oscillatory instability, and 
that liquid peroxide explosions may be initiated at an oscillatory instability 
rather than by classical thermal ignition. Similar results have been obtained 
using parameters for the decomposition of cumene hydroperoxide in equations 
(3-4) (Ball, 2010). 

Is it realistic to model a reacting volume inside a storage tank — or, for that 
matter, a peroxide bomb — as a well-stirred flow reactor? Yes, on a timescale 
over which the reacting volume remains relatively constant and gradientless 
relative to the much faster rate of reaction. For the purposes of this analysis in 
which the focus is on the dynamics we can circumscribe a reacting volume in 
which the spatial gradients are insignificant in comparison to the time evolution, 
and therefore can be neglected. In this case the CSTR paradigm is appropriate. 
If this approximation does not hold, then we are free to reduce the circumscribed 
volume until it does. There is nothing particularly artificial or manipulative in 
doing this; it is just a simplest case scenario for which the powerful tools of 
stability and bifurcation theory can be applied. Much of the heat transfer 
would be convective rather than conductive, and on convective timescales the 
approximation does not hold — but that is for a separate study. 

5.1. Nanoscale aspects of oscillatory thermal instability 

Since this is a book about all things nano, it is pertinent to elucidate the 

mechanism of oscillatory thermal instability on the nanoscale. But before we 
zoom in to nano spatial scales we need to discuss the characteristics of relaxation 
oscillators in general and examine the two time scales of the relaxation oscillation 
solutions of equations 3 and 4. 

Relaxation oscillators are often and easily implemented in electrical and 
electronic circuits, but any dissipative dynamical system with nonequilibrium 
steady states and two or more dynamical state variables has the potential to 
exhibit relaxation oscillations. In general they are not smoothly sinusoidal in 
form. Instead, relaxation oscillators are characterised by energy dissipation and 
energy accumulation occurring on different timescales. There may be relatively 
slow dissipation until the system reaches a threshold state at which the internal 
energy rapidly and nonlinearly increases, or rapid dissipation followed by slow 
but accelerating increase in internal energy. This latter two-time dynamics 
occurs with the thermochemical oscillator in equations 3 and 4. 

The two time scales in each period evident in figure 5 are shown in figure 9. 
On the 'fast' time scale, left, during the few seconds between the temperature 
maximum and about 990 s the concentration remains close to zero. On the 'slow' 
time scale, right, the system barely heats at all over the long interval, although 
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Figure 9: Left: a 'fast' or dissipative time interval. Right: a 'slow' time interval. 



evidently some reaction takes place because reactant accumulates nonlinearly, 
reaches a maximum, then declines slowly before the exponential spike. 

We can extract these two time scales approximately from equations 3 and 
4. Define 'fast' or 'stretched' or dissipative time as r' = r/e and recast the 
equations with t' as the independent variable: 

^=s{-xe-'/- + f{l^x)) (5) 

^^:ce-^/- + ef{^Ua~u)~£{u~Ua). (6) 

For s sufhciently small we have dx/dr' ~ and x ~ ~ since the reactant 
is almost fully depleted, and du/dr' ~ —£{u — Ua), so that u decays as u w 
exp(— £r')(uo — Ua) + Ua = exp(— £r')(uinax ~ Ua) + Ua since the temperature 
amplitude is at a maximum when a; = at r^^ax = ''o = 0- As r' becomes 
'large' u — > Ua', in figure 9 left this occurs around ~1000s. However reactant is 
gradually accumulating and the system evolves in 'slow' time. 

Over much of the 'slow' time interval the system behaviour can be under- 
stood by dividing equation 4 through by e. At low temperatures and with large 
e we then have du/dr « — /(w — Ua), which gives u ~ exp(— /t)(uo — Ua) + Ua, 
and since uq = Ua on this time scale we have u Ua- In figure 9 this approx- 
imation holds over much of the 'slow' time interval: the internal energy of the 
system increases but the temperature remains almost constant, although slowly 
increasing. During this time reactant accumulates, but reaction does take place: 
X evolves as x sa /(I — exp(— 6r)/6, where b = exp(— I/uq). 

If reaction takes place in 'slow' time but the system is barely heating up, 
where does the heat of reaction go? Now we focus on the nanoscale. 

The answer is that the heat is stored in the internal motions of the reactant 
mixture molecules. Referring to the definition of the dimensionless group e in 
table 3 we see that e is large if the specific heat capacity Cr of the reaction 
mixture is large. The specific heat capacity of a substance is a function of the 
number of degrees of freedom of motion that are available to its constituent 
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particles. A monoatomic perfect gas has only the kinetic energy of each atom 
so translational motion in three dimensions is the only motion it can undergo. 
The equipartition theorem tells us that the three translational contributions to 
the constant volume molar heat capacity C'^J = (3/2)i? = 12.47 J/(molK), thus 
the kinetic energy contribution does manifest as temperature change. 

Polyatomic molecules have many additional degrees of freedom though, and 
physical properties of liquids arc governed much more by the potential energy of 
the system than the kinetic energy. Potential energy is stored in intramolecular 
rotational and vibrational degrees of freedom and in the vibrational intermolec- 
ular force potential. The mean translational energy of each liquid molecule is 
the same as that for gases, (3/2) fcT, where k is Boltzmann's constant, and since 
Cp = Cy + R we can take the molecular translational contribution to the molar 
heat capacity of a liquid as 

C*' = 20.78Jmor^K-^ 

Since the molar heat capacity of MIC is 111.76 J mol~^K~^ and that of water 
is 75.29 JmoP^K^^, on a molar basis the MIC-H2O liquid system can store a 
large amount of heat in the internal rotational modes and in the vibrational 
intermolecular potentials. (The internal vibrational modes are not excited at 
this 'slow' temperature.) 

The specific heat capacity of MIC, as opposed to the molar heat capacity, 
is not particularly large, but that of water is anomalously high due to strong 
intermolecular forces. For the reaction mixture it is it is large enough to depress 
the first and third terms on the right hand side of equation 4 relative to the 
second term at low temperature. The rapid temperature spike and transition 
to 'fast' time occurs because the activation energy is high enough to make the 
reaction very temperature-sensitive. In other words, the Arrhenius term in 
equation 4 can take over very rapidly once the reaction zone reaches a certain 
temperature. At the temperature peak the system enters the 'fast' dissipative 
time interval, since the reactant is fully depleted. The specific heat capacity 
cannot affect the dynamics on this time scale. 

This analysis suggests that containment systems with large heat capacity 
for thermoreactive liquids may not suppress the oscillatory instability, although 
they would certainly lengthen the 'slow' timescale. This may or may not be 
a good thing. A lengthy 'slow' interval may provide enough time before the 
thermal runaway to deactivate or quench the system in other ways. On the 
other hand, it might give a false sense of security: 'Nothing has happened over 
the last four hours, we might as well leave it. ' 

6. Conclusions 

1. The CSTR paradigm was applied to investigate the thermal stability of 
MIC hydrolysis and the thermal decomposition of triacetone triperoxide 
in solution. 
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2. Stability analyses of the steady state solutions of the dynamical model 
found that in both cases thermal runaway occurs due to the hard onset of 
a thermal oscillation at a subcritical Hopf bifurcation. Classical thermal 
ignition at a steady state turning point does not occur in these systems 
and over the thermal regime of interest they are dominated by oscillatory 
instability. 

3. This non-classical oscillatory thermal misbehaviour may be generic in liq- 
uid thermoreactive systems where the specific heat capacity and activation 
energy are high. These results provide new information about the cause 
of thermal runaway that may inform improved designs of storage systems 
for thermally unstable liquids and better management of organic peroxide 
based explosives. 

Appendix 

The following brief account of the Bhopal disaster has been compiled from 
the following sources: Forman (1985); Weir (1987); Shrivastava (1987); Lep- 
kowski (1994) and Abbasi and Abbasi (2005). 

The Union Carbide plant at Bhopal carried out the production of carbaryl, 
an agricultural insecticide that has been used widely throughout the world since 
1945. Methyl isocyanate, a low-boiling point, highly reactive and extremely 
toxic liquid used in the synthesis of carbaryl, was stored in an underground 
stainless steel tank (Tank 610) which was encased in a concrete shell. The 
temperature of the 41 tonnes of MIC in Tank 610 was 12-14° C rather than the 
recommended 0-4° C because the refrigeration unit had been non-operational 
for several months. On the evening of December 2 1984 a worker had been sent 
to hose out a nearby tank. The hose was left running unattended, and it is 
believed that a faulty valve allowed entry of water into the connected Tank 610. 
(Union Carbide disputes this, asserts that nothing was wrong with its equipment 
and procedures, and argues that sabotage by a disaffected employee must have 
led to the disaster.) By 11:30 pm, when workers detected lachrymose whiffs of 
leaking MIC, water had been running into Tank 610 for several hours. Although 
a slow rise in temperature and pressure in the tank had been noted, the early 
signs of trouble were not acted upon. Shortly after 11:30 pm the contents of 
the tank reached thermal criticality and began escaping as vapor from the flare 
tower. 

Downwind of the flare tower lay the crowded suburbs and shantytowns. Most 
of the fluid in the tank streamed from the tower then drifted low over the city 
and sank and seeped in deathly mist in lungs in eyes, a period to sleep and 
swift arrest of retreat. About 4000 lives were claimed immediately, many tens 
of thousands through the subsequent days and months and years lost their lives 
or their health to the poison's efi'ects, and the dead are still being counted. 

Acknowledgement: The author is recipient of Australian Research Council Fu- 
ture Fellowship FT0991007. 
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